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Role  of  Spin  Momentum  Current  in  Magnetic  Non-Local  Damping  of 

Ultrathin  Film  Structures 

G.  Woltersdorf,  R.  Urban,  and  B.  Heinrich 
Simon  Fraser  University,  8888  University  Dr.,  Burnaby,  BC,  V5A  1S6,  Canada 

Non-local  damping  was  investigated  by  Ferromagnetic  Resonance  (FMR)  using  ultrathin  magnetic 
single  and  double  layer  structures  prepared  by  Molecular  Beam  Epitaxy  (MBE).  The  double  layer 
structures  show  magnetic  damping  which  is  caused  by  spin  transport  across  a normal  metal  spacer 
(N).  In  double  layer  structures  a thin  Fe  layer,  FI,  was  separated  from  a thick  Fe  layer,  F2,  by  a 
Au(001)  spacer.  The  interface  magnetic  anisotropies  separated  the  FMR  fields  of  FI  and  F2  by  a big 
margin  allowing  one  to  investigate  FMR  in  FI  while  F2  had  a negligible  angle  of  precession,  and  vice 
versa.  The  Fe  films  in  magnetic  double  layers  acquire  non-local  interface  Gilbert  damping.  It  will 
be  shown  that  the  precessing  magnetic  moments  act  as  spin  pumps  and  spin  sinks.  This  concept 
was  tested  by  investigating  the  FMR  linewidth  around  an  accidental  crossover  of  the  resonance 
fields  for  the  layers  FI  and  F2.  There  is  another  possible  mechanism  for  non-local  damping  which  is 
based  on  a ’’breathing  Fermi  surface”  of  the  spacer.  The  temperature  dependence  of  the  non-local 
damping  indicates  that  this  mechanism  is  weak  in  Au  spacers.  Surprisingly  the  Au  spacer  acts  as  an 
additional  impedance  for  the  spin  pump  mechanism.  Finally,  it  will  be  shown  that  electron-electron 
correlations  in  a Pd  spacer  can  lead  to  a significant  enhancement  of  the  non-local  damping. 


I.  INTRODUCTION 


The  small  lateral  dimensions  of  spintronics  devices  and  high  density  memory  bits  require  the  use 
of  magnetic  metallic  ultrathin  film  structures  where  the  magnetic  moments  across  the  film  thickness 
are  locked  together  by  the  intra  layer  exchange  coupling.  Since  spintronics  and  high  density  magnetic 
recording  employ  fast  magnetization  reversal  processes  it  is  important  to  understand  the  spin  dynamics 
and  magnetic  relaxation  processes  of  multilayers  in  the  nano-second  time  regime.  The  spin  dynamics  is 
described  by  the  Landau  Lifshitz  Gilbert  (L.L.G.)  equation  of  motion  / 
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where  7 is  the  absolute  value  of  the  electron  gyromagnetic  ratio,  Ma  is  the  saturation  magnetization  and 
G is  the  Gilbert  damping  parameter.  The  effective  field  H e//  is  given  by  the  derivatives  of  the  Gibbs 
energy,  U,  with  respect  to  the  components  (Mx,  My,Mz)  of  the  magnetization  vector  M(f)  [1],  The 
second  term  in  eq.  1 represents  the  well  known  Gilbert  damping  torque. 


II.  NON-LOCAL  DAMPING:  EXPERIMENT 


The  role  of  non-local  damping  was  investigated  in  high  quality  crystalline  Au/Fe/Au/Fe(001)  struc- 
tures grown  on  GaAs(OOl)  substrates,  see  details  in  [2-5].  In-plane  Ferromagnetic  Resonance  (FMR) 
experiments  were  carried  out  using  10,  24,  36  and  72  GHz  systems  [4,  6]. 

Single  Fe  ultrathin  films  with  thicknesses  of  8,  11,  16,  21,  and  31  monolayers  (ML)  were  grown  directly 
on  GaAs(OOl)  and  covered  by  a 20  ML  thick  Au(001)  cap  layer  for  protection  in  ambient  conditions.  FMR 
measurements  were  used  to  determine  the  in-plane  four-fold  and  uniaxial  magnetic  anisotropies,  K\  and 
Ku,  and  the  effective  demagnetizing  field  perpendicular  to  the  film  surface,  4irMeff,  as  a function  of  the 
film  thickness  d [6].  The  magnetic  anisotropies  were  well  described  by  the  bulk  and  interface  magnetic 
properties,  respectively  [2].  The  thin  Fe  films  which  were  studied  in  the  single  layer  structures  were  grown 
again  as  a part  of  magnetic  double  layer  structures.  The  thin  Fe  film  (FI)  was  separated  from  the  second 
thick  Fe  layer  (F2)  by  a Au(001)  spacer  (N)  of  a variable  thickness  between  8 to  100  ML.  The  magnetic 
double  layers  were  covered  by  a 20  ML  Au(001)  capping  layer.  The  thickness  of  the  Au  spacer  layer  was 
always  smaller  than  the  electron  mean  free  path  (38  nm)  [7],  and  hence  allowed  ballistic  spin  transport 
between  the  magnetic  layers. 


67 


The  interface  magnetic  anisotropies  separated  the  FMR  fields  of  FI  and  F2  by  a big  margin  (~1 
kOe),  see  Fig.l  allowing  one  to  carry  out  FMR  measurements  in  FI  with  F2  possessing  a small  angle 
of  precession  compared  to  that  in  FI,  and  vice  versa.  The  thin  Fe  film  in  the  single  and  double  layer 
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FIG.  1:  (a)  The  FMR  fields  at  24  GHz  in  the  layer  Fl(16Fe,  shown  by  (o))  and  layer  F2(40Fe,  (*))  in 
20Au/40Fe/16Au/16Fe/GaAs(001)  (integers  represent  the  number  of  atomic  layers)  as  a function  of  the  angle  ip 
between  the  applied  field  and  the  in-plane  [100]  crystallographic  direction.  A large  in-plane  uniaxial  anisotropy 
field  in  FI  leads  to  an  accidental  crossover  at  =115  and  150  Deg.  Notice  that  the  FMR  fields  get  locked  together 
by  the  spin  pumping  effect  at  the  accidental  crossover.  Away  from  the  crossover  the  resonance  fields  are  separated 
by  as  much  as  ten  FMR  linewidths. 

(b)  The  FMR  linewidths  corresponding  to  FI  in  Fig.  1(a).  The  measured  and  calculated  FMR  signals  were  an- 
alyzed using  two  Lorenzian  lineshapes.  The  Lorenzian  peaks  were  characterized  by  their  amplitudes,  resonance 
fields  and  linewidths.  The  solid  lines  were  obtained  from  calculations  using  eqs  7.  Note  that  the  FMR  linewidth 
for  FI  first  increases  before  it  reaches  its  minimum  value  corresponding  to  that  of  a single  layer  structure. 


structures  had  the  same  FMR  field  showing  that  the  static  interlayer  exchange  coupling  in  the  double 
layer  [2]  through  the  Au  spacer  was  negligible.  The  absence  of  static  interlayer  coupling  is  due  to  short 
atomic  terraces  [1]. 

The  results  are  as  follows:  (a)  The  FMR  linewidth  in  the  thin  films  FI  always  increased  in  the  presence 
of  F2,  and  vice  versa;  (b)  The  additional  FMR  linewidth,  A Hadd,  followed  an  inverse  dependence  on 
the  thin  film  thickness  d\,  see  [2];  and  (c)  the  additional  FMR  linewidth  A Hadd  in  both  the  parallel  (H 
in-plane)  and  perpendicular  (H  perpendicular  to  the  plane)  FMR  configuration  was  linearly  dependent 
on  the  microwave  frequency  with  no  zero  frequency  offset  [5]  . The  additional  Gilbert  damping  for  the 
16  ML  thick  Fe  film  was  found  to  be  only  weakly  dependent  on  the  crystallographic  direction,  with  the 
average  value  Gad<i  — 1-2  x 10s  s_1.  This  is  comparable  to  the  intrinsic  Gilbert  damping  in  the  single  Fe 
film,  Gint  =1.4xl08  s-1. 


III.  THEORY  OF  NON-LOCAL  DAMPING 


Tserkovnyak  et  al.  [8]  showed  that  the  interface  damping  can  be  generated  by  spin  current  from  a 
ferromagnet  (F)  into  the  adjacent  normal  metal  reservoirs  (NM).  The  spin  current  is  generated  by  a 
precessing  magnetic  moment.  A precessing  magnetization  at  the  F/NM  interface  acts  as  a ’’peristaltic 
spin  -pump”.  The  direction  of  the  spin  current  is  perpendicular  to  the  F/NM  interface  and  points  away 
from  the  interface  towards  NM.  The  spin  momentum  which  is  carried  away  by  the  spin  current  is 
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FIG.  2:  A cartoon  representing  the  dynamic  coupling  between  two  magnetic  layers  which  are  separated  by  a 
non-magnetic  spacer  N.  (a)  represents  two  magnetic  layers  with  different  FMR  fields.  FI  is  at  resonance,  and  F2 
is  nearly  stationary.  A large  grey  arrow  in  the  normal  spacer  describes  the  direction  of  the  spin  current.  The 
dashed  lines  represent  the  instantaneous  direction  of  the  spin  momentum.  For  small  angle  of  precession  they 
are  nearly  parallel  to  the  transverse  rf  magnetization  component  shown  in  short  solid  arrows.  FI  acts  as  a spin 
pump,  F2  acts  as  a spin  sink,  (b)  represents  a situation  when  FI  and  F2  resonate  at  the  same  field.  Both  layers 
act  as  spin  pumps  and  spin  sinks.  In  this  case  the  net  spin  momentum  transfer  across  each  interface  is  zero.  No 
additional  damping  is  present. 


where  m is  the  a unit  vector  in  direction  of  M.  The  spin  current  can  cause  magnetic  damping.  Ar  is  the 
interface  scattering  parameter,  and  for  F films  thicker  than  spin  coherence  length,  7t /(kp  — kp),  given  by 

Ar  = \Yl\r1mn-rinn\2  , (3) 

m,n 

where  rj^n  are  the  electron  reflection  matrix  elements  at  the  NM/F  interface  for  the  spin  up  and  down 
electrons,  k are  the  Fermi  k vectors  for  majority  and  minority  electrons,  respectively.  The  sum  in  Ar 
is  close  to  the  number  of  the  transverse  channels  in  NM  [9].  The  sum  is  given  by 

T = i=°-85’*w'  <4> 


where  S is  the  area  of  the  interface,  kF  is  the  Fermi  wavevector  and  n is  the  density  of  electrons  per  spin 
in  NM  [9].  Brataas  et  al.  [9,  TO]  showed  that  Ar  can  be  evaluated  from  the  interface  mixing  conductance 
G-\i  [11].  Ar  = -^G][  = Sg^i,  where  represents  ’’dimensionless  interface  mixing  conductivity”. 

The  generated  spin  current  propagates  though  the  normal  metal  spacer  N,  and  is  deposited  at  the 
N/F2  interface.  It  was  shown  [10,  12]  that  the  transverse  component  of  the  spin  current  in  N is  entirely 
absorbed  at  the  N/F2  interface,  see  Fig.2a.  For  small  precessional  angles  the  spin  current  is  almost 
entirely  transverse.  This  means  that  the  N/F2  interface  acts  as  an  ideal  spin  sink,  and  provides  an 
effective  spin  brake  for  the  precessing  magnetic  moment  in  FI  [4,  5, 13].  The.  spin  momentum  }spin  which 
is  carried  away  by  the  spin  current  has  the  form  of  Gilbert  damping.  The  resulting  Gilbert  damping  is 
given  by  the  conservation  of  the  total  spin  momentum 
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where  Mj0t  is  the  total  magnetic  moment  in  FI.  After  simple  algebraical  steps  one  obtains  an  expression 
for  the  dimensionless  spin  pump  contribution  to  the  damping,  aap, 
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where  d\  is  the  thickness  of  FI,  is  the  dimensionless  mixing  conductivity,  and  Gsp  is  the  spin  pump 
Gilbert  parameter,  g is  the  electron  g-factor.  The  inverse  dependence  of  asp  on  the  film  thickness  clearly 
testifies  to  its  interfacial  origin.  The  layers  FI  and  F2  act  as  mutual  spin  pumps  and  spin  sinks.  For 
small  precessional  angles  the  eq.  of  motion  for  FI  can  be  written  as  [4] 
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where  Mi  is  the  magnetization  vector  of  FI,  mi;  are  the  unit  vectors  along  Mi;,  and  d%  is  the  thicknesses 
of  FI.  The  exchange  of  spin  currents  is  a symmetric  concept  and  the  eq.  of  motion  for  the  layer  F2  is 
obtained  by  interchanging  the  indices  1 2.  The  third  and  fourth  terms  in  eq.  7 represent  the  spin 

pump  and  spin  sink  of  FI.  The  fourth  term  is  generated  by  the  spin  pump  from  F2.  For  clarity  it  is 
worthwhile  to  point  out  that  the  signs  (+)  and  (— ) in  the  third  and  fourth  terms  in  eq.  7 represent  the 
spin  current  directions  (FI  — > F2)  and  (F2  — > FI),  respectively. 

The  spin  pump  model  is  a rather  exotic  theory  to  those  who  are  working  in  magnetism.  One  would 
expect  that  there  is  a direct  connection  to  a more  common  concept  which  is  applicable  to  magnetic  mul- 
tilayers. The  obvious  choice  is  interlayer  exchange  coupling.  In  the  past  interlayer  exchange  interaction 
was  treated  only  in  the  static  limit  [14].  Recently  it  was  shown  that  its  dynamic  part  can  create  magnetic 
damping  [4,  13,  15].  One  can  show  that  the  spin  pumping  theory  is  equivalent  to  the  dynamic  response 
of  the  interlayer  exchange  coupling.  Since  the  dynamic  coupling  is  derived  from  the  rf  susceptibility  one 
can  account  for  electron-electron  correlation  effects  in  the  normal  metal  spacer  [13].  It  has  been  shown 
[15]  that  the  Gilbert  damping  is  enhanced  by  the  square  of  the  Stoner  factor  Se  — [1  — UN(Ep)]^1, 


G?_hd  = Gs-dS%,  (8) 

where  U is  the  screened  interatomic  Coulomb  interaction  and  N(Ep ) the  electron  density  of  states,  per 
atom,  at  the  Fermi  level  in  NM. 

It  is  worthwhile  to  realize  that  the  s — d exchange  relaxation  mechanism  also  applies  to  bulk  ferro- 
magnets,  and  was  evaluated  by  Heinrich  et  al.  [16,  17].  The  Gilbert  damping  in  this  case  is  given 
by 

GMk  = f (9) 


where  xp  is  the  Pauli  susceptibility  and  r»/  is  the  spin  flip  relaxation  time  of  itinerant  electrons  in  the 
ferromagnet.  It  should  be  noted  that  1 /rsf  in  metals  is  proportional  to  the  square  of  the  spin  orbit 
interaction  [16,  17].  Using  xp  from  Kriesman  and  Callen  [18]  and  rs/  from  the  spin  diffusion  length  in 
Current  Perpendicular  to  Plane  (CPP)  GMR  measurements  one  obtains  for  the  bulk  Gilbert  damping 
G = 5xl06s-1  and  1 xl08s-1  for  Co  and  permalloy  (Py),  respectively,  see  the  details  in  [17].  This 
contribution  is  small  in  Co  but  it  explains  the  intrinsic  damping  in  Py.  Fe  is  expected  to  behave  like  Co. 

One  expects  that  there  has  to  be  an  additional  mechanism  which  depends  explicitly  on  rs/.  Its  origin 
can  be  understood  in  the  following  way.  The  mechanism  of  the  interlayer  exchange  coupling  lies  in 
the  itinerant  nature  of  the  electron  carriers  [13].  The  energy  of  electrons  depends  on  the  instantaneous 
orientation  of  the  magnetic  moments,  and  consequently  the  occupation  number  n k)<r  of  electronic  states 
having  energy  e k,<r  changes  during  precession  of  the  magnetization,  and  this  results  a ’’breathing  Fermi 
surface”,  this  concept  was  also  used  in  bulk  materials  [19,  20].  However,  this  redistribution  cannot 
be  achieved  instantaneously.  The  time  lag  between  the  instantaneous  exchange  field  and  the  induced 
moment  in  the  spacer  is  described  by  the  transverse  spin  relaxation  time,  ra/,  which  is  proportional  to 
the  momentum  relaxation  time  entering  conductivity  [21].  In  the  limit  of  slow  precessional  motion  the 
effective  damping  field  can  be  evaluated  [4,  13], 
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where  the  sum  is  carried  out  per  unit  area.  In  this  case  the  damping  field  is  proportional  to  the  spin 
relaxation  time  rsf.  This  mechanism  is  explicitly  dependent  on  the  conductivity  and  it  represents  a 
different  contribution  to  the  non-local  damping  compared  to  the  spin  pumping  mechanism  which  is 
independent  on  r,/. 


IV.  DISCUSSION  OF  THE  RESULTS 

The  spin  pumping  and  breathing  Fermi  surface  theories  predict  Gilbert  damping,  which  results  in  a 
strictly  linear  dependence  of  A Hadd  on  the  microwave  frequency.  This  is  indeed  experimentally  observed 
over  a wide  range  of  microwave  frequencies  [3,  4]. 

The  validity  of  the  spin-pumping  theory  can  be  tested  by  comparing  calculations  using  eqs  7 with  the 
experimental  results.  Fig.2  shows  two  extreme  situations.  In  Fig.  2(a)  the  FMR  fields  in  FI  and  F2  are 
separated  by  a big  margin.  In  Fig.2(b)  the  FMR  fields  are  same.  In  (a)  one  expects  the  full  contribution 
from  the  non-local  damping.  AHadd  for  FI  and  F2  should  scale  with  their  respective  1/d  terms.  In  (b)  the 
situation  is  symmetric;  the  net  spin  momentum  current  trough  both  interfaces  is  zero,  and  the  additional 
damping  is  absent.  This  behavior  is  well  demonstrated  in  Fig.lb.  The  good  agreement  between  theory 
and  experiment  clearly  shows  the  validity  of  the  spin  pumping  and  spin  sink  concept  which  is  described  by 
eqs  7.  Even  in  the  absence  of  static  interlayer  exchange  coupling  the  magnetic  layers  are  coupled  by  the 
dynamic  part  of  the  interlayer  exchange.  The  spin  sink  effect  at  the  N/F  interface  starts  to  be  inefficient 
only  when  the  thickness  of  the  normal  metal  spacer  N becomes  comparable  to  the  spin  diffusion  length 
[22] . Since  the  resistivity  ratio  in  our  structures  is  a factor  of  two,  see  Fig.4  the  spin  diffusion  length  in 
Au  is  at  least  of  the  order  of  100  nm.  The  static  interlayer  exchange  coupling  vanishes  in  our  samples  due 
to  interface  roughness  on  a length  scale  of  a mere  10  ML  (2  nm).  One  should  point  out  that  when  the  N 
metal  spacer  thickness  starts  to  be  comparable  to  the  spin  diffusion  length  then  the  N spacer  absorbs  a 
part  of  the  spin  current  [23,  24]. 
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FIG.  3:  The  additional  FMR  linewidth,  AHadi,  as  function  of  the . Au  spacer  thickness  in  20Au/8- 

c_  27Au/16Fe/GaAs(001).  The  samples  with  the  Au  spacer  thickness  between  8 to  16ML  were  grown  separatelly. 
The  samples  with  the  spacer  thickness  between  17  and  27ML  were  grown  using  a wedged  Au(001)  layer.  One  can 
see  some  degree  of  variation  in  A Hadd  with  a period  of  ~ 3ML.  These  oscillations  are  marginal. 

The  quantitative  comparison  with  predictions  of  the  spin  pumping  theory  is  very  favorable.  First 
principles  electron  band  calculations  [11]  resulted  in  g jj  ~l.llxl015cm-2  for  a clean  Cu/Co(lll)  interface. 
By  scaling  this  value  to  Au  using  eq.4  one  obtains  Gsp  =1.37xl08s_1  which  is  only  14%  higher  than  the 
value  measured  at  room  temperature.  This  is  surprising  agreement  considering  the  fact  that  calculations 
of  the  intrinsic  damping  in  bulk  metals  have  been  carried  out  over  the  last  three  decades,  and  yet  they 
have  not  been  able  to  produce  a comparable  agreement  with  experiments  [17]. 

The  breathing  Fermi  surface  contribution  to  the  Gilbert  damping  is  proportional  to  the  electron  relax- 
ation time  Tsf  of  the  N metal  spacer,  see  eq.10.  Since  this  contribution  is  using  the  concept  of  interlayer 
exchange  coupling  one  can  expect  some  degree  of  oscillatory  behavior  with  a changing  spacer  thickness. 
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FIG.  4:  The  temperature  dependence  of  the  additional  FMR  linewidth,  A Hadd  for  3 differ- 
ent samples:  (a)  20Au/40Fe/14Au/16Fe/GaAs(001),  (★)  20Au/40Fe/40Au/16Fe/GaAs(001),  and  (.) 

20  Au/40Fe/100Au/16Fe/GaAs(001) . The  temperature  dependence  of  the  sheet  conductivity,  cr,  for  the  (•)- 
sample  is  shown  in  the  dashed  line.  Note  that  the  temperature  dependence  of  A Hadd  is  linear  and  weak  for  all 
samples.  One  can  fit  the  temperature  dependence  by  eq.  11  . 


FIG.  5:  The  slope,  (dA Hadd/dT),  of  the  additional  FMR  linewidth,  as  function  of  the  spacer  thickness.  The 
slope  was  obtained  for  the  samples  shown  in  Fig.4.  The  solid  line  is  a spline  fit  to  guide  the  reader’s  eye. 


The  measurements  shown  in  Fig.  3 indicate  that  this  behavior  is  only  weakly  present.  A better  test  of 
the  breathing  Fermi  surface  contribution  can  be  carried  out  by  measuring  the  temperature  dependence 
of  the  non-local  damping.  One  expects  proportionality  with  the  sheet  conductance  (rs/  ~ Torb  ~ a ) of 
the  N spacer.  The  temperature  dependence  of  the  sheet  conductance  and  the  additional  FMR  linewidth 
is  shown  in  Fig.  4.  The  temperature  dependence  od  A Badd  is  linear  and  obviously  does  not  scale  with 
the  conductivity.  This  clearly  indicates  that  the  strength  of  the  breathing  Fermi  surface  contribution  is 
unimportant  in  the  Fe/Au/Fe(001)  samples.  Note,  that  the  slope  of  A Hadd  as  a function  of  temperature 
increases  with  an  increasing  spacer  thickness.  One  can  fit  the  temperature  dependence  by  the  linear 
function 

A-ffadd  = Slope  • T + AHadd(T  = 0 K).  * (11) 

Note,  that  AHadd(T  = 0 K)  is  equal  for  all  measured  Au  spacers.  Fig.  5 shows  the  fitted  slope  as  function 
of  the  spacer  thickness.  One  can  see  that  the  slope  saturates  for  higher  spacer  thicknesses.  This  implies 
that  a part  of  the  Au  spacers  at  finite  temperatures  decreases  effectiveness  of  the  spin  pump  mechanism 
and  acts  as  an  additional  resistance  in  series  with  the  interface  mixing  resistance  l/t/n.  This  decreases 
the  spin  current  reaching  the  N/F2  interface.  The  additional  resistance  in  Au  decreases  with  decreasing 
temperature  and  becomes  zero  for  T = 0 K.  Since  the  spin  diffusion  length  in  Au  is  much  bigger  than 
the  Au  spacer  thickness  no  additional  resistance  in  Au  is  expected.  The  linear  temperature  dependence 
of  AiJatjc;(T)  can  be  indicative  that  the  Au  resistance  follows  the  temperature  dependence  of  inverse 
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FIG.  6:  (a)  The  FMR  fields  of  the  two  Fe  layers  in  the  GaAs/16Fe/9Pd/40Au/40Fe/20Au(001)  sam- 
ple. (•)  and  (★)  correspond  to  16Fe  and  40Fe,  respectively.  The  measurements  where  carried  out  at 
24GHz.  (b)  shows  the  FMR  linewidth  for  16Fe  in  single  [GaAs/16Fe/9Pd/20Au(001)]  (o),  and  double  layer 
[GaAs/16Fe/9Pd/40Au/40Fe/20Au(001)]  (•)  sample.  We  did  not  analyze  the  FMR  data  around  the  accidental 
crossover  of  resonance  fields  at  110  and  165  degrees,  respectively.  Around  these  points  the  two  FMR  peaks  resulted 
in  a broad  and  asymmetric  spectrum  not  allowing  one  to  carry  out  a reliable  analysis. 


susceptibility  {1/x  ~ T)  of  ” loose  spins”  [1,  25].  The  impurity  Fe  atoms  (loose  spins)  can  be  interdiffused 
up  to  certain  depth  (see  saturation  of  this  effect  with  increasing  Au  spacer  thickness  in  Fig.  3)  around  the 
Fl/N  interface  by  surfactant  Arsenic  during  the  deposition  of  Au  on  FI  [7].  Using  this  picture  one  can 
argue  that  the  temperature  dependence  of  Af?a(jrf(T)  in  Fig.  3 could  be  a consequence  of  ’’loose  spins” 
which  contribute  to  the  impedance  of  the  spin  pumping  mechanism.  It  is  worthwhile  to  note  that  the 
measured  value  of  Gsp  — 1.4xl08s-1  at  T = 0 K is  very  close  to  that  estimated  from  the  spin  pumping 
theory,  Gap  =1.37xlCrs-1. 

The  dynamic  exchange  coupling  theory  [13],  includes  the  Stoner  enhancement  factor,  see  eq.8.  Our 
recent  results  using  20Au/4Pd/[Fe/Pd]5/14Fe/GaAs(001)  single  and  20Au/40Fe/40Au/4Pd/[Fe/Pd]s- 
/14Fe/GaAs(001)  double  layer  samples,  see  [13],  show  a strong  evidence  for  the  Stoner  enhancement 
factor.  This  structure  incorporates  a magnetic  [Fe/Pd]s  superlattice  with  5 repetitions  of  a [IFe/lPd]  unit 
cell.  The  N metal  spacer  is  4Pd40Au(001).  The  additional  FMR  linewidth  along  the  cubic  crystallographic 
axes  was  enhanced  by  as  much  as  a factor  of  four,  see  [4,  13]  and  is  significantly  bigger  than  that  expected 
from  the  simple  spin  pumping  mechanism.  It  turns  out  that  one  can  get  the  same  enhancement  in  non- 
local damping  using  a simple  GaAs/16Fe/9Pd/40Au/40Fe/20Au(001)  double  layer  structure,  see  Fig.  6 
where  the  [Fe/Pd]s  superlattice  is  replaced  by  5 atomic  layers  of  Pd.  In  this  case  the  single  layer  structure 
is  GaAs/16Fe/9Pd/20Au(001)  and  shows  no  enhanced  damping,  see  Fig.  6.  The  metallic  Pd  is  known  to 
exhibit  a strong  Stoner  enhancement  in  the  dc  susceptibility  [26] . The  Fermi  surface  of  Pd  is  not  spherical 
and  can  lead  to  a significant  angular  dependence  of  the  Stoner  enhancement  factor  which  satisfies  the 
(001)  symmetry  of  our  samples,  see  Fig.  6.  These  results  clearly  show  that  electron  correlation  effects  in 
normal  metal  spacers  have  to  be  considered. 


V.  CONCLUSIONS 

We  have  shown  that  non-local  damping  by  the  transfer  of  spin  momentum  can  be  realized  in  magnetic 
multilayer  films.  This  effect  is  significant  in  ultrathin  films.  Theoretical  models  were  presented  for  the 
non-local  damping.  It  has  been  demonstrated  that  the  non-local  interface  Gilbert  damping  in  magnetic 
multilayers  is  well  described  by  the  concept  of  spin  pumps  and  spin  sinks.  It  has  been  shown  that  this 
effect  is  directly  related  to  the  dynamics  of  the  interlayer  exchange  coupling.  By  proper  engineering  of 
multilayer  structures  one  can  create  magnetic  damping  which  significantly  surpasses  that  in  the  bulk 
materials. 
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